
restore,'./irdc_dist_model/beta_cube.sav',/ver

myps,'./irdc_dist_model/analysis_plots/beta_plot.eps',xsize=10,ysize=5.0,/cmyk
multiplot_xm,[2,1],xgap=0.00,/xmargin

med_pch = median(pchoo_arr,dim=1)
plot,med_pch,ratios,/ylog,xr=[0.5,0.999],/xst,thick=3,$
     ytit='Probability Ratio of Peaks',$
     xtit='Integrated Probability on Chosen Side of d!dtan!n'
vline,0.683,/ylog
vline,0.955,/ylog
vline,0.997,/ylog

;; vline,/h,2,thick=5,color=cgColor('BLK5')

cgLoadct,13,/silent
cgColorbar,/right,/vert,range=[0.0,1.0],$
           position=[0.92, 0.16, 0.95, 0.90],charsize=1.0,$
           title='Cumulative Fraction as a Function of Ratio'

multiplot
cgLoadct,13,/silent
plotimage,rhist_arr,xst=4,yst=4,range=[0.0,1.0],$
          min_dpi=100

cgContour,rhist_arr,levels=[0.15,0.50],/over,c_thick=3,$
          c_annot=['15%','Median']

cgContour,pchoo_arr,levels=[0.6827,0.9545,0.9973],/over,c_thick=3,$
          c_color='White',charsize=0.8,$
          c_annot=['68.27%','95.45%','99.73%']

axis,xaxis=0,xr=minmax(alpha),/save,/xst,xtit=cgSymbol('beta')
axis,xaxis=1,xtickformat='blank_axis'
axis,yaxis=0,yr=minmax(ratios),/ylog,/save,/yst,$
     ytit='';'Probability Ratio of Peaks'
axis,yaxis=1,ytickformat='blank_axis'

;; vline,/h,2,thick=5,color=cgColor('BLK5')
;; vline,2,/ylog,thick=5,color=cgColor('BLK5')


myps,/done,/mp
end
